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Abstract 

A model with a singular forward scattering amplitude for particles with 
opposite spins in d spatial dimensions is proposed and solved by using the 
bosonization transformation. This interacting potential leads to the spin- 
charge separation. Thermal properties at low temperature for this Luttinger 
liquid are discussed. Also, the explicit form of the single-electron Green func- 
tion is found; it has square-root branch cut. New fermion field operators are 
defined; they describe holons and spinons as the elementary excitations. Their 
single particle Green functions possess pseudoparticle properties. Using these 
operators the spin-charge separated Hamiltonian for an ideal gases of holons 
and spinons is derived and reflects an inverse (fermionization) transformation. 
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It was suggested |]] that the properties of normal state of high-temperature supercon- 
ductors are properly described by Luttinger liquid, where the spin and charge degrees of 
freedom are separated. In one-dimensional systems this phenomenon is well understood [0 . 
However, in two and three dimensions the present understanding of the spin-charge sepa- 
ration is rather poor. In this letter we formulate and solve exactly a d-dimensional model 
exhibiting the spin-charge separation, as well as discuss its thermal and dynamic properties. 

A natural approach to study spin-charge decoupling phenomena is the bosonization trans- 
formation, generalized recently to the multidimensional space situation ||. Here we adopt 
the operator version of the bosonization developed in Ref. @]. The starting assumption in 
this method is the existence of the Fermi surface (FS) defined as a collection of points at 
which the momentum distribution function has singularities at zero temperature (T = 0). 
These points are parameterized by vectors S and T, which label a finite and a locally flat 
(rectangular in shape) mesh of grid points on FS with spacing A <C kp between them 
Introducing coarse-grained density fluctuation operators J CT (S,q), defined in boxes centered 
at each FS point and having surface area A d_1 and the thicknesses A/2 both above and be- 
low it, one can transform the effective Hamiltonian for interacting fermions into an effective 
Hamiltonian for free bosons. Explicitly, it takes the general form 

H = \ E E E r -'( S > T, q)J„(S, q)J^(T, -q), (1) 

1 S,T q crcr' 

where r CT0 .'(S, T, q) = VF(S)^5 a y5^ + -^V^ CT /(S, T, q) is the positive defined matrix el- 
ement. The first term corresponds to the kinetic energy part of the original fermionic 
Hamiltonian with linearized dispersion relation close to FS, whereas the second term is the 
effective (low-energy) interaction between the particles with spins a and a'. The geometrical 
factor fl = A d ~ 1 (-^) d depends on the system dimension d. 

The explicit expression for V aa >(S,T,q) is generally derived by transforming out the 
high energy modes in the fermionic Hamiltonian ||. Obviously, this procedure can also 
change the Fermi velocity vf{S). Therefore, we take vf{S) as an effective value obtained 
after removing the high-energy degrees of freedom. As shown below, we can characterize 

2 



the universal properties of fermions knowing only the asymptotic behavior of the interaction 
potential V CT(7 /(S, T, q) in the thermodynamic limit. 

If the system is invariant under the time reversal, the interaction part must be explic- 
itly symmetric under this operation, which means that V^ CT /(S,T, q) = S, — T, — q). 
Furthermore, if FS is also invariant under the reflections S — > — S etc., the last condition 
becomes V^ CT /(S, T, q) = V^'(S,T,q). In that case V aa '(S, T, q) depends only on the rel- 
ative orientation of the spins a and a'; there are only two independent components: V arj 
for parallel spins and V a „ for antiparallel spins. It is convenient to introduce the symmetric 
and the antisymmetric combinations: \^ C,S (S, T, q) = ^{V aa (S, T, q) ± V^(S, T, q)), where 
c and s superscripts corresponds to "±" signs, respectively. Correspondingly, we define the 
currents J CjS (S, q) = J](S, q) ± Jj,(S, q)), which describe the charge and the spin density 
fluctuations, respectively. Then, the original Hamiltonian (p]) takes the form 

H = E ^MSlk^ql^-ql + kEHSJ.ql^qlJ^T.-q). 

a=c,s z s q S,T q 

(2) 

The a = c term describes the dynamics of the charge density fluctuations in the system, 
whereas the a = s term deals with the longitudinal spin density fluctuations. There are no 
terms which mix the degrees of freedom (i.e. ~ J c ■ J s ) because the Hamiltonian is assumed 
to be invariant under the spin flip (i.e. J c — > J c , J s — > — J s ). We can check out that for 
the noninteracting case, i.e. for V a = 0, both the spin and the charge density fluctuations 
propagate with the same velocity vp(S). The commutation relation for density fluctuation 
operators take the following form [J a (S, q), Jp(T, p)] = ^«/3^s/r^p+ q ,o ^ 1 ' ^s, where a, (3 = 
c, s. Thus, the two branches of density fluctuations are independent of each other. The 
commutation relations become equivalent to those obeyed by the bosonic harmonic-oscillator 
creation and annihilation operators after rescaling them by the factor on the right hand side, 
i.e. by defining the creation (a+) and the annihilation (a a ) operators according to 

J Q (S,q) = 6(n s ■ q)\/tth s ■ q a Q (S, q) + 6(-n s ■ q)yj-£lfi s ■ qa+(S,-q), (3) 
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where 0(x) is the step function. Then a a (S, q), dp(T, p) = Sg^S^ 5 a p. I 11 terms of the 
bosonic harmonic-oscillator creation and annihilation operators the Hamiltonian (|2]) is 



H= E E^E {^s • q)0(fi T • q) ^(vf(S) ■ q)(v F (T) • q) x (4) 



Oc=C,S s,T q 




We have shown previously |||7] that a universal behavior of the system takes place and 
depends on how the interaction part of the effective Hamiltonian (|j) behaves in the scaling 
(thermodynamic) limit A — ► 0. If the interaction in this limit has a singular power-law 
behavior, i.e. V a ^(S, T, q) scales for S — ► T as A d ~ 1 /A ri , then: (i) for r\ < d — 1 the Landau 
Fermi liquid (FL) fixed point is stable; (ii) for rj > d — 1 the statistical spin liquid (SSL) 
is the stable fixed point; (iii) the case rj = d — 1 leads to the Luttinger liquid (LL) type 
of behavior. Properties of FL in d = 3 and 2 are well known SSL was discussed in 
Refs. Here we concentrate on the properties of the spin-charge separated LL in an 

arbitrary spatial dimension. 

To examine the basic properties of Luttinger liquid we assume the following form of the 
effective interaction in ([!]): 



y TT (S,T,q) = / TT (S,T,q), 



(5) 



and 
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where / CT(T /(S, T, q) and g(S,q) are nonsingular functions. In this manner, we model the 
situation with a divergent forward scattering amplitude (taking place for S = T and q = 0), 
which leads to the LL fixed point. Also, since we examine only the low-energy limit, we can 
expand the nonsingular part according to: g(S, q) = go(S) + <?(S)ns ■ q + and omit the 
higher-order terms. Substituting this expansion into (|B|) we obtain in the thermodynamic 
limit (A — > 0) the Hamiltonian with two branches of free bosons excitations, each with a 
different form of the kinetic energy. Namely, defining Vp S (S) = Vp(S)±g(S), the Hamiltonian 
(|) simplifies to 

H = E E K(S) • q) a+(S, q) a a (S, q). (7) 

a=c,s S,q>0 

The charge and the spin fluctuations propagate in the system with different velocities and 
express the separation of the corresponding degrees of freedom. One of the velocities dimin- 
ishes and the other increases. In the extreme case one of them vanishes transforming into a 
soft mode, signalling a phase transition in the corresponding channel, c or s. The character 
of this transition will be determined by the regular part. The spin-charge decoupling takes 
always place in the one- dimensional systems of interacting fermions in the low energy limit 
0. In a system of higher dimension the potential must be singular. Obviously, this is a 
simplified model. More appealing form would be Vj^(S,T, q) ~ |s-T|*+|q| 7 ' > which clearly 
diverges in the forward direction (for q — > 0). Such a singular effective potential in d = 2 was 



discussed in Refs. |]10|-|12| in the context of the high-temperature superconductors and the 
fractional quantum Hall effect. Our choice (§) is modeled by the most divergent term of this 
potential, and physically means that the fermions with antiparalel spins interact through 
the forward scattering processes along the radial FS direction only. 

Since the Hamiltonian (0) is diagonal we can calculate the internal energy of the system 
and then the specific heat. The energy is 

^ = g.+ EEE e ;!g ) q _ 1 , (8) 

a=c,s S 1 

where we utilize the fact that now we are dealing with bosons. Eq is the ground state 
energy of the initial noninteracting system. This term must be incorporated because 



in the bosonization procedure the energy is measured with respect to the Fermi energy. 
Also, the chemical potential \x does not appear in (§) because the number of bosons 
is not conserved. In other words, those fields describe the system particle-hole excita- 
tions. This is also the reason why those bosons cannot condense. For an isotropic sys- 
tem the velocities do not depend on the index labeling the point of FS, i.e. f^(S) = v F , 
and we find that U = ^(k B T) 2 (J^J p^+i) ^f" 1 W + ^q- Hence, the specific heat 
is Cy = ja^f = ^(A;b) 2 T (^r) d f^TTy^" 1 E« ^r- Introducing the density of states for 
the charge and the spin excitations on FS: p a (ep) = (j^j r(rf/2+i) ^f" 1 ^) we have that 
Cy = j(^B) 2 rEa Pa(tF) = IllT. The spin-charge separated liquid has a linear spe- 
cific heat, with jll ~ 55- + ^r- The linear specific heat is thus a general characteristic 
following from the existence of the FS, independently of the statistical properties of the 
particles. In the limiting case vp = v s F = vp we recover the FL result. Similarly, we 
find the free energy in the form F = E — ^(/c^T) 2 J2 a Pa{ e F), an d the low-T entropy 
S = = ^(^j^Eft Po(^f), which coincides with the specific heat. Thus the low-T 
thermal properties of the present Luttinger spin-charge separated liquid are very similar to 
those of free fermions. The only difference is in the form of the density of states at FS. 
However, the dynamic properties in the LL case are quite unique, as we discuss next. 
We define the fermion correlation function as 

G>(S, x, t > 0) =< V.(S, x, t)^+(S, 0, 0) > . (9) 

To derive an explicit form of this function we substitute the bosonized Fermi field op- 
erators H Vv(S,x) = y / J e iks - x e i ^^ (s ' x) 6(S), and utilize the identity e A e B =: e A+B : 
e < AB ~hi A2 + B2 )> ) where :: means that : e A+B : is a normal ordered product of operators. 
Then, we have that G>(S,x, t) = ^-exp (^^G B (S : x, t)J , where G B is expressed via the 
Bose fields, namely G B (S,x,t) =< (0 c (S,x,t) + a<p s (S, x, t))(0 c (S, 0, 0) + ct0 s (S,O,O)) > 
— < (0 C (S,O,O) + cr0 s (S, 0, 0)) 2 > . Next, using the Heisenberg representation for the 
boson field operators <f> a (S,x,t), and subsequently decomposing them into the Fourier 
components, we find the explicit form of the boson correlation function: G B (S,x.,t) = 
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~ In ( ^ ns x "f t+ "^" s x v F t+ia) \ ^ jj ence ^ fermion correlation function is 

4-7T V \ la ) } 



e iks-x 



g>(s, x, t) = in (io) 

'fis ■ x — f^t + iaJhs ■ x — v F t + ia 



It is independent of the spin index a. We see that instead of a quasiparticle pole, taking 
place in the FL case, we have now a branch cut ranging from v s F to v c F . This branch cut 
survives when we transform G^(S, x, t) into G^(S, k, a;) as can be easily seen by decomposing 
the argument into normal and transverse parts: k • x = k(hs • x + t • x), where tk is the 
transverse part of k, and noting that the part hs ■ x can be transformed in the same manner, 



as in the d — 1 case [13]. The t ■ x part has a trivial form. This means that the analytic 
character of the LL Green function is the same in both d — 1 and d > 1 cases. This universal 
character follows from the relation G > (S, k, uj) = G > (S, n s • k, u;)5 t .k,o- This very significant 
result tells us also that there are no quasiparticle excitations having a direct relation to 
the noninteracting particles in this charge-spin separated system. In other words, when we 
put a single electron forming a wave packet on FS, it dissociates into many wave packets 
propagating with velocities ranging form v s F to v c F . This nonperturbative result means that 
there is no one-to-one correspondence between the dynamics of LL liquid system and the 
system of non-interacting fermions. Note that the distribution function is in our model 
situation at T = a step function with ci jump cit as in the FL case. The inclusion 
the nonsingular part of the interaction will change this step distribution for a finite- volume 
system. 

The fundamental question arises if we can still define a proper fermionic pseudoparticles 
in this spin-charge separated liquid. To construct such a state we consider the field operators 
ip c and ip s , defined through the following fermionization transformation 

lMS,x) = ^e lks -V^-( s ' x )0(S). (11) 

The operators ip a (S, x) obey proper anticommutation relations. By the procedure similar to 
that employed in deriving (0), we now have the correlation function in the new fermionic 
variables in the form 



G>(S, x, t > 0) =< ^ Q (S, x, t)iP+(S, 0, o) >= in- — e ' ks " ■ , , (12) 

ns ■ x — Vpt + la 

which has simple poles. Also, one can show that [^(S, x), J dyJ a (S,y)] = t/>+(S,x), which 
means that ^>+(S, x) changes the total number of either charge (a = c) or spin (a = s) of 
the system at the FS point S by one unit. Hence, the operators ([□]) represent new fermionic 
pseudoparticles for the interacting non-Fermi liquid. Those single-particle excitations are 
called the holons for a = c and the spinons for a — s. They are the only single-particle 
excitations across FS in our model system. Since they obey the fermion statistics, the 
contribution to the specific heat is of the fermionic type. In other words, the spin-charge 
separated liquid is composed of ideal gases of spinons and holons. They represent the exact 
eigenstates of the system. Effectively, we have the following Hamiltonian for noninteracting 
fermion pseudoparticles 

H= E 5>S(S) /rfx^(S,x)(^Ws,x). (13) 

a=c,s S \ * / 

Since we have treated the interaction between fermions exactly in the thermodynamic limit 
(A — > 0), the spectrum of our new noninteracting Hamiltonian ( Jl3| ) is exactly the same 
as that of the interacting Hamiltonian ([!]). To conclude, the two branches of elementary 
excitations in this LL can be represented either as boson (collective) or as fermion (single- 
particle) excitations. 

Having determined the excitation spectrum of the singular part of interaction, which 
resulted in the spin-charge separated liquid, we can now include in (|13D the nonsingular 
parts of the interaction, transformed to the fermionic representation (|TTD . Since for the 
former part we determined an exact state of the system, we can now treat the residual 
interaction among them as a perturbation, i.e. regard the system of interacting holons and 
spinons as being in one-to-one correspondence to the system of noninteracting holons and 
spinons. This statement can be proved by referring to the Gell-Mann and Low theorem 
which means that the evolution operator is well defined in all orders, since the interactions 
among the holons and the spinons are nonsingular functions. In other words, the adiabatic 
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"switching on" procedure is justified. Here we discuss only a semiclassical approach to the 
interacting holons and spinons in the bosonic language. The equations of motion for J a (S, q) 
operators are 

J a (S, q, t) = v£(S) • q J a (S, q,t) + q-n s A^ 1 (^) £ /«( S " T ) J «( S ' * *)> ( 14 ) 
where we supposed that the interaction f a does not depend on q. In the semiclassical 
approach we take the expectation value of J a , i.e. define u a (S, q, t) =< J a (S, q, t) >. 
This quantity measures the shape deformation of FS with respect to the ground state form. 
Furthermore, focusing our attention on the single Fourier mode u a (t) = e~ luJt u a , we find 
the integral equation for the collective excitation spectrum amplitude u a = u a (Q, q), in the 
form 

(v pq cos 9 — u)u a (fl,q) = qvp cos 9 / -—F a (Cl,Cl')u a (Cl,',q), (15) 

where Sd = J dfl, F a (fl,Q') = p(eF)f a {&,Q'), and Q is the solid angle. This is an equation 
of motion for either the holon (a = c) or the spinon (a = s) sound-wave amplitudes. The 
introduced bosons represent the sound waves propagating around FS, here characterized 
by uj and q. Eq. fllBD gives both the stable solution for collective modes and the solution 
with the imaginary frequency. Solution of Eq. fll5|) is analogical to that considered in the FL 
theory ||; it will not be discussed in detail here. 

In summary, we presented a model with the spin-charge separation in the space of arbi- 
trary dimensions, as well as have discussed some of its basic properties. The next step would 
require a careful analysis of the nonsingular part of the interaction for finite- volume systems. 
In particular, the most important question is to construct a theory of interacting holons and 
spinons at low energies (in the spirit of the Landau FL), including also the effects coming 
from the presence of an applied magnetic field. One should also examine the stability of this 
liquid against charge or spin-density wave formation. 
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